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ABSTRACT 

We model the radio, optical, and X-ray emission for the afterglows of GRB 980703, 990123, 990510, and 
991216, within the framework of relativistic jets, to determine their physical parameters. The models that yield 
acceptable fits to the data have jet energies mostly between 10 50 to 10 51 erg and initial opening angles between 1° 
and 4°. The external medium density is uncertain by at least one order of magnitude in each case, being around 
10~ 3 cm' 3 for GRB 980703 and 990123, - 10" 1 cm" 3 for GRB 990510, and - 3 cm" 3 for GRB 991216. If 
the jets are uniform (i.e. there are no angular gradients of the energy per solid angle) then the 20 keV - 1 MeV 
radiative efficiency during the GRB phase must have been at least 2-3% for GRB 990510, 20% for GRB 990123, 
and 30% for GRB 991216. 

Subject headings: gamma-rays: bursts - ISM: jets and outflows - methods: numerical - radiation mechanisms: 
non-thermal - shock waves 



1. INTRODUCTION 

There are two basic quantities one needs to try understand- 
ing the nature of any astronomical source - the distance and the 
energy of the source. For the long duration GRBs, lasting more 
than 10 seconds, the former is well established to be cosmolog- 
ical. However, the energy associated with the GRBs remains 
uncertain. 

The efficiency of producing 7-ray emission in the generally 
accepted internal shock model (Meszaros & Rees 1994), which 
can explain the observed temporal variability, is less than a 
few percent (Kumar 1999, Lazzati, Ghisellini & Celotti 1999, 
Panaitescu, Spada & Meszaros 1999, see also Beloborodov 
2000). This suggests that the total energy in the explosion is 
larger than the energy observed in 7-rays by a factor of 20-50. 

The goal of this work is to infer the physical parameters of 
the ejecta, including their energy and the external medium den- 
sity, from modeling of radio, optical and AT -ray data for GRB 
afterglows with known redshifts. The modeling is carried out in 
the framework of collimated ejecta interacting with an isotropic 
external medium. The model is described in <|2] and the results 
of the numerical calculations for individual afterglows are given 



2. DESCRIPTION OF THE MODEL 

In calculating the jet dynamics, we assume that the energy 
and baryon density within the ejecta do not have an angular 
dependence, and that the external medium is isotropic. We 
also assume that, at any time, the physical parameters and bulk 
Lorentz factor V are the same in the entire swept up external 
gas. We include the effect of radiative energy losses on the jet 
dynamics. 

For the calculation of synchrotron emission, we assume a 
tangled magnetic field and that the electrons accelerated at 
shock have a power-law distribution. The electron distribu- 
tion resulting from the continuous injection at shock and adi- 
abatic+radiative cooling is approximated by broken power-law, 
with a break at the minimum random Lorentz factor of the 



freshly injected electrons and another one at the cooling elec- 
tron Lorentz factor. In calculating the received flux, the swept 
up gas is approximated as a surface, i.e. the thickness of the 
emitting shell is ignored. The effect of this surface curvature 
on the photon arrival time and energy is taken into account. 

2.1. Dynamics 

The interaction between the relativistic ejecta that generated 
the GRB with the external gas continuously decelerates the jet 
and heats the newly swept up gas. Assuming that the heated gas 
has a uniform temperature, equal to that of the freshly shocked 
fluid, the total energy of the GRB remnant is: 



E{r) = m(r)(T 2 - 1) + m {T - 1) 



(1) 



where too is the mass of the ejecta, whose Lorentz factor 
and energy at the end of the GRB phase are Tq and Eq = 
(To — 1)tooc 2 , respectively. 

In the above equation, to is the mass of the swept up external 
gas, given by 

dm(r) = S}(r) p(r)r 2 dr , (2) 

where p(r) cx r~ s is the external medium density (s = for 
homogeneous gas, s = 2 for a wind ejected at constant speed 
before the release of the ultra-relativistic ejecta), and 0(r) is 
solid angle of the jet. 

The jet half opening angle increases due to the lateral spread- 
ing of the jet at the local sound speed c s : 

rd9 = c s dt' = r- 1 c s d£ /Q6 , (3) 

t' and tiai, being the time measured in the ejecta comoving and 
laboratory frames, respectively. The speed of sound is (see 
Huang et al. 2000) 



7(7 -iK 2 7(7-i)(r-i) 



p' + ju' 



i + 7 (r-i) 



(4) 



where v! = (T — l)p' and p' are the comoving internal energy 
and rest mass densities, respectively, and 7 is the adiabatic in- 
dex. In the relativistic limit 7 = 4/3 and c s — c/\/3, while for 
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non-relativistic speeds 7 = 5/3 and c s = (v5/3)u, where v is 
the radial expansion speed. We relate the adiabatic index with 
T through a simple formula, which has the above asymptotic 
limits. 

The energy losses through synchrotron and inverse Compton 
emissions is given by 



AE 

dtlab 



a e c m(r) 
67r m rj 



"I'M 



B 2 (Y + 1) 7 W( 7 ), (5) 



where B is the magnetic field intensity, Y the Co mpt on param- 
eter, M is the normalized electron distribution (§ |2.2| ), and 7 is 
the electron random Lorentz factor. 

Equations ([[]), (||), (|^) and (||) are solved numerically, subject 
to the boundary conditions: T(0) = T , m(0) = 0, 9(0) = 9o, 
and E(0) = E . 

2.2. Electron Distribution and Spectral Breaks 

The magnetic field intensity is parameterized relative to its 
equipartition value 



B 2 = 8vrp'c 2 (r- l)e B = 32tteb p{r)c 2 (T - 1 



7T + 1 

'FT' 



(6) 

where p 1 is the comoving frame rest-mass density. 

The distribution of the electrons accelerated by the forward 
shock and injected in down-stream is assumed to be a power- 
law of index p 



M(7) °t 7 p , 7i < 7 < 7m , 



(7) 



where 7i is the minimum, injected electron Lorentz factor, pa- 
rameterized relative to its value at equipartition, 



7* 



(8) 



and 7m is an upper limit, determined by the conditions that the 
acceleration timescale of such electrons does not exceed the ra- 
diative losses timescale, and that the total energy in the injected 
electrons does not exceed a certain fraction e of the available 
internal energy. The former condition leads to 



3e 



1 



g(T e B(Y + 1) 



1/2 



(9) 



where n g is the ratio of the acceleration timescale to the gyra- 
tion time. The latter condition can be written as 



m , 



7dM(7)=en*p(T-l) , 



(10) 



where Mi is normalized (to unity), and m e and m„ are the elec- 
tron and proton mass, respectively. Equation (|0|) leads to an 
algebraic equation which can be solved numerically. The upper 

limit 7m is the minimum between 7^ and yff above. Unless 
rig is larger than about 10 3 , the upper limit given in equation 

([)]) is sufficiently high that the synchrotron emission from 7^ - 
electrons is above the soft X-ray domain. However if p < 2.5 
and e is not much larger than e e , the upper limit given by equa- 
tion ( |To| ) may be sufficiently low to yield a break of the after- 
glow emission at the frequencies of interest (X-rays and even 
optical). For numerics we shall use n g = 10 and e = 0.5, 



the latter corresponding to equipartition between electron and 
protons. 

The distribution of cooled electrons is a power-law of index 2 
if the electrons are cooling faster than the dynamical timescale, 
and a power-law steeper by unity than the injected distribution 
in the opposite case (Sari, Piran & Narayan 1998). Therefore 
the electron distribution resulting from injection at shock and 
radiative cooling is 



7V(7) oc 



7" 



7c < 7 < Tt 
7» < 7 < 7m 



(11) 



for fast cooling electrons ( 7c < 7 0, and 



7V( 7 ) oc 



7 



-v 

-(p+i) 



7i < 7 < 7c 
7c < 7 < 7 M 



(12) 



for slow cooling electrons ( 7 j < j c ). In equations ( pT[ ) and 
(p"2|), 7c is the cooling electron Lorentz factor, defined by the 
equality of its radiative cooling timescale with the dynamical 
timescale: 



7c = 67T- 



m e c 



1 



(T e t'B 2 {Y + l) 



(13) 



The Compton parameter Y is given by 
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7 2 dA/"( 7 ) 



(14) 



where 7 m = min( 7 i, 7c ) and r e is the optical thickness to elec 
tron scattering: 

ct e m(r) 



m p ll(r)r 2 



(15) 



The Klein-Nishina effect reduces the inverse Compton losses 
above an electron Lorentz factor •jkn approximated as the ge- 
ometric mean of i) the electron Lorentz factor for which scat- 
tering of the synchrotron photons emitted by such an electron 
occurs at the Klein-Nishina limit and of ii) the electron Lorentz 
factor for which scattering of the synchrotron photons emitted 
by 7 m -electrons is at the same limit. The comoving frame syn- 
chrotron characteristic frequency for an electron of Lorentz fac- 
tor 7 is 



A-r) = 



16 m e c 



BY 



(16) 



We take into account the Klein-Nishina reduction by calculat- 
ing the integral in equation ( |i4| ) up to jkn if Ikn < 1m, and 
by switching off the inverse Compton losses above jkn in the 
integral given in equation (||). 

The synchrotron self-absorption frequency in the fluid rest 
frame is at 1/ = v'(j a ) with 7 a given by (see Panaitescu & 
Kumar 2000) 



5e r e 



3/10 



-1/2 



(17) 



This equation is valid only if 7a < 7r) 
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2.3. Received Flux 

The synchrotron spectrum is approximated as piece-wise 
power-law (see Sari et al. 1998) with breaks at the injection, 
cooling, and absorption breaks given by equations (fig), ([$]), 
©.and©. 

To calculate the afterglow flux seen by the observer, we con- 
sider that the emitting shell is infinitely thin and that the ob- 
server is located on the jet axis. Consider an annular region of 
area 5 A — 2tt r 2 Sp, with p = cos u, where u> is the polar angle, 
measured relative to the jet axis. The energy emitted in the co- 
moving frame per unit time and frequency is SL' V , = P' v , H e SA, 
where P' v , is the radiative comoving power per electron and 
E e is the electron surface density. The infinitesimal comoving 
energy emitted per solid angle (\ / A-K)8L' ll ,dv l dt 1 is relativisti- 
cally beamed toward the observer by a factor T> 2 and boosted 
in frequency by a factor T>, where T> = [T(l — /3/i)] -1 , with 
= v/c. Therefore the infinitesimal flux dF v received by the 
observer at frequency v = T>v' , during At = c~ 1 r5p, satisfies 



dF„ dv St 



1 + z 



V 3 PLT, e dt' dv' 5 A , 



(18) 



where P' v , is the sum of synchrotron and inverse Compton emis- 
sions, z is the afterglow redshift, and is the luminosity dis- 
tance. We assume a Universe with Hq = 65 kms _1 Mpc _1 , 
9o = 0.1, and A = 0. 

The flux received by the observer at time t is that given by 
equation dl8f), integrated over the entire evolution of the source. 
Using dt' = di; a h/r = dr/(/3cT) and relating the electron sur- 
face density to the jet mass and area, S e = m(r)/ (m p fiV 2 ), 
equation (|8j) leads to 



2 m p d? L 



dr 



P' v , (r)m(r) 



/3r 3 (l-/3^) 2 rO(r) 



(19) 



with /i given by the condition that light emitted from location 
(r, p) arrives at observer at time t = ti a i — (r/c)/i. Thus equa- 
tion ( 19 ) takes into account the spread in the photon arrival time 
due to the spherical curvature of the jet surface. In all our cal- 
culations it is assumed that the observer is on the jet axis. 

3. ANALYTICAL CONSIDERATIONS 

So far there are five afterglows (990123, 990510, 991216, 
000301c, 000926) for which a break in the optical emission has 
been identified. In all these cases the break is seen at or after 
t = 1 day. Within the framework of uniform ejecta interacting 
with isotropic media, there are two possible causes for such a 
break: the passage of a break frequency (injection, Vi, cooling, 
v c , or that due the upper cut-off of the electron distribution, 
vm), or the edge of the jet becoming visible to the observer 
(plus the changing dynamics due to the lateral spreading of the 
jet). 

One can show that, within a factor of order unity, the break 
frequency vi is the same for both types of external medium: 



10 13 {z + lf'%% 



54 c e, 



1/2 .-3/2 tt 
-l £ B,-2 t d nz 



(20) 



where £q — 47r(£'o/rio) is the jet isotropic equivalent energy, 
td is observer time in days, and the usual notation A n = l0~ n A 
was used. Equation ( p0[ ) shows that, unless the isotropic equiv- 
alent energy exceeds 10 56 erg and the magnetic field is close to 
equipartition (e e cannot be much higher than 0.1, as the frac- 
tional energy in electrons must be below unity), vi is below the 



optical range at t <; 1 day. Therefore it is very unlikely that the 
optical light-curve breaks are due to the passage of z/j through 
the observational band. Moreover, if this were the case, then, at 
times before the light-curve break, the temporal index a of the 
light-curve decay, F v [t) oc t~ a , would be at most 1/4, which 
is much smaller than the observed a's. 

3.1. Passage of the Cooling Break 

We consider here the afterglow emission at early times, when 
the effects due to collimation of ejecta are negligible, but suf- 
ficiently large that v. L < v. In this case, the afterglow light- 
curves for slow cooling electrons (z/j < v c ) are given by (see 
Panaitescu & Kumar 2000) 



F v>Vi oc 



t -(3p-2)/4 V C <V,Y <\ 

r (3p/4)+l/(4-p) Vc < ^ y > 1, 2 <p < 3 
t -(3p/4) + l v c <v, Y >1, 3<p 

(21) 

for a homogeneous external medium (s = 0), where the last 
two rows represent the case when the electron cooling is domi- 
nated by inverse Compton scatterings. For a wind-like medium 
(s = 2) and slow cooling electrons 



t -(3p-l)/4 
i-(3p-2)/4 
f -(3p/4)+p/(8- 
f -(3p/4)+3/2 



2p) 



Y < 1 

Y > 1, 2<p<3 

Y > 1, 3<p 

(22) 

The second row of equations ( j2Tj ) and ( 22 ) also gives the light- 



v < v c 

Vc < V, 
V c < V, 



curve for Vi < v and fast cooling electrons. 

The temporal evolution of the cooling break frequency v c for 
s = and slow cooling electrons is given by 



Y < 1 

,/ -v <; i -(8-3p)/(8-2p) Y > 1 , 2 < p < 3 

t 1 / 2 Y > 1 , 3 <p 

For s = 2 and slow cooling electrons 



(23) 



t i/ 2 

f(3p-4)/(8-2p) 
£5/2 



Y < 1 

Y > 1 

Y > 1 



2 <p < 3 

3 < p 



(24) 



The first row also gives the evolution of v c for fast cooling elec- 
trons, in which case Y is time-independent. 

For a homogeneous medium, equation ([0]) shows that v c in- 
creases in time if the electron cooling is dominated by inverse 
Compton and if p > 8 /3. From equation (pl|), the passage of 
v c through the observational band changes the light-curve de- 
cay index by 



(H 



1/4 Y < 1 

(8-3p)/(16-4p) y>l,2<p<8/3 

(3p-8)/(16-4p) y>l,8/3<p<3 

1/4 Y > 1 , 3 <p 



(25) 

Note (Aa) c > 0, i.e. the passage of v c always steepens 
the light-curve decay, even if v c increases in time, and that 
(Aa) c < 1/4. In the case where v c is above optical and be- 
low X-ray, the temporal indices of the X-ray and optical light- 
curves differ by a x — a Q = (Aa) c > if u c decreases in time, 
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and by a x — a = — (Aa) c < if v c increases in time. Thus, 
for s — 0, the X-ray emission decays faster than the optical one 
if Y < lorF > land p< 8/3. 

For a wind-like medium equation ( p4| ) shows that v c always 
increases in time. From equation the passage of v c yields 



for s = Q and 



r 1/4 

(A«) c = (3p-4)/(16-4p) 
I 5/4 



y < i 

Y> 1 , 2 <p < 3 

y > i , 3 < P 



(26) 

Note that 1/4 < (Aa) c < 5/4. For t/ Q < v c < v x , the X-ray 
and optical indices differ by a x — a = — (Aa) c < 0. Hence, 
for s = 2, the X-ray emission always decays slower than the 
optical one. 

3.2. Collimation ofEjecta 

If the ejecta is collimated, the decay of the afterglow emis- 
sion steepens around the time tj when T6 = 1, due to the al- 
tered jet dynamics and that the observer sees the edge of the jet. 
For s = 



tj = 1.2 (z + 1) (£ o ,540g,_ino 1 ) V3 da y 



(27) 



The coefficient above was determined numerically from the ar- 
rival time of the photons moving toward the observer along the 
jet axis. Photons emitted from other regions on the jet surface 
arrive later by a factor up to ~ 4. 

Around tj the jet dynamics changes from a quasi-spherical 
expansion with T cx r~ 3 / 2 oc t~ 3 / 8 for s — (T cx r -1 / 2 oc 
t~ x l A for s = 2) to a sideways expansion characterized by 
r oc e~ kr cx £ -1 / 2 (Pvhoads 1999). During the lateral spreading 
phase {t > tj) the cooling frequency evolution is 



t l 



Y < 1 



V c CX I <2(p-2)/(4-p) Y >l . 2 < V < 3 



V 

y > 1 , 3 <p 



(28) 



assuming slow cooling electrons (y. L < v c ). Then it can be 
shown that, at t > tj, the light-curve is given by 



t~ p 

£-p+(p-2)/(4-p) 

t -CP-i) 



any v , 
v c < v 

< V 



v, 



Y < 1 

y > i 
y > i 



2 <p < 3 

3 <p 

(29) 

Evidently, as the source slows down, the Y parameter eventu- 
ally falls below unity and the last two cases given in equation 
( P9[ ) approach asymptotically F v cx t~ p . The results given in 
equations ( |28| ) and (g9|) ignore multiplying terms that are pow- 
ers of the jet radius r, which increases logarithmically with the 
observer time. With the same approximation, they also hold for 
a jet interacting with a pre-ejected wind. Furthermore, these 
results are accurate only at times when the afterglow is very 
relativistic. From numerical calculations we found that, for the 



first case given in equation ( 29 ), the decay index a is approxi- 
mated by p with an error less than 10% if T 10. 

Using equations (21 ), (^2|), and (^) it can be shown that, for 
v > Vi, the magnitude of the break due to collimation of ejecta 



(A«)j 



(p + 3)/4 
(p+2)/4 
p/4 + l/(4-p) 
(p + 4)/4 



v <v c 

v c < v, Y < 1 

u c < v, Y > 1, 2 < p < 3 ' 
i/ c < i/, Y > 1, 3 < p 

(30) 



(Aa)j 



(p+l)/4 i/ < 

(p + 2)/4 v <v,Y<\ 
p/4 + p/(8-2p) i/ c < v, Y > 1, 2 <p < 3 
(p + 6)/4 i^c <v, Y > 1, 3 <p 

(31) 

for s = 2. The finite opening of the ejecta yields Aq = 3/4 for 
s = and Aq = 1/2 for s = 2 (Panaitescu, Meszaros & Rees 
1998) when the jet edge becomes visible, the remainder of the 
steepening being due to the sideways expansion of the jet (Ku- 
mar & Panaitescu 2000). Equations ((3^) and (|T]) show that, if 
p > 2, (Aa)j > 1 for s = and (Aa) 3 > 3/4 for s = 2 . 

3.3. What Can We Infer from the X-ray 
and Optical Decay Indices ? 

The most important difference between a break due to pas- 
sage of v c and one due to collimation of ejecta is the chro- 
maticity of the former and the achromaticity of the latter. This 
would be the best criterion to distinguish between them if opti- 
cal and X-ray observations spanning the same 1-2 decades in 
time, around the time when the break i s se en, are available. 

The analytical results presented in §3.1 and fjfO allow us to 



draw some conclusions even when the existence of simultane- 
ous X-ray and optical light-curve breaks cannot be clearly es- 
tablished. Equations (25 ), (pq), (|30|), and (pl|) show that optical 



break magnitudes Aq < 3 /4 can be produced only by the pas- 
sage of the cooling break, while Aa > 5/4 can be due only 
to collimation of ejecta. The caveat of this criterion is that, as 
shown by Kumar & Panaitescu (2000), for collimated ejecta, 
the completion of most of (Aa)j is spread over at least one 
decade in observer time for s = and over at least two decades 
for s = 2. Therefore observations spanning a shorter time range 
will underestimate the true magnitude of the jet edge break, par- 
ticularly in the case when the observer is not located close to the 
jet axis. 

The results presented in §3.1 and §3.2 also lead to the fol- 



lowing criteria for determining the location of v c relative to the 
optical and X-ray domains, and for identifying the type of ex- 
ternal medium, from the optical and X-ray decay at the same 
time: 

1. if a a — a x , then either v c is not between the optical and 
X-ray domains, or it is in this range but is quasi-constant 
in time. The latter may occur either in collimated ejecta 
during the sideways expansion phase, or before jet edge 
effects are important, in s = models, if the inverse 
Compton cooling is the dominant process and if p ~ 8/3 
(see eq. [||). 

2. if a x > a (i.e. the X-ray decay is faster than in the 
optical), then v < v c < v x and the external medium 
is homogeneous. Furthermore either the electron cool- 
ing is synchrotron dominated or, if Compton-dominated, 
p < 8/3. 

3. if a x < a a , then v a < v c < v x and either the exter- 
nal medium is a wind-like or it is homogeneous and the 
electron cooling is inverse Compton-dominated (Y > 1), 
withp > 8/3. 

Using the observed spectral slope in the optical range and the 
relative intensity of the X-ray and optical emission, one can 
eliminate some of the above cases, and further reduce the num- 
ber of potentially good models for a given afterglow. 
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4. INDIVIDUAL AFTERGLOWS 

The modeling of GRB afterglow light-curves is carried out 
by solving numerically equations (|]), (g|), (||), and ( B), t o deter- 
mine the dynamics of the afterglow, and equation (|l9|) to cal- 
culate the observed flux. The six unknown parameters £o, 9q, 
n, p, e e , and eb are determined by minimizing the x 2 between 
observed and model fluxes at the frequencies where most of the 
data is available. The jet initial Lorentz factor Tq does not af- 
fect the afterglow emission if its evolution is quasi-adiabatic, 
because it drops out from the calculation of the jet Lorentz fac- 
tor r as a function of observer time. When radiative losses 
are taken into account, Tq may have an effect on the after- 
glow emission, as the magnitude of these losses depends on Tq. 
However, the effect is weak and cannot be used to set significant 
constraints on Tq. For simplicity, we keep a fixed To — 500 in 
all numerical calculations. 

In this work we restrict our attention to four afterglows 
for which radio, optical, X-ray light-curves and redshifts are 
available: GRB 980703, GRB 990123, GRB 990510, and 
GRB 991216, leaving out GRB 970508, whose optical light- 
curve cannot be entirely explained within the framework of our 
model, as it exhibited a brightening after 1 day, indicating a 
possible delayed energy injection, or fluctuations of the energy 
release parameters e e , eb or of the external density n. 

The near-infrared (NIR) and optical magnitudes are con- 
verted to fluxes using the characteristics of photometric bands 
published by Campins, Rieke, & Lebofsky (1985), and 
Fukugita, Shimasaku & Ichikawa (1995). We assumed a 5% 
uncertainty in the magnitude-flux conversion. The NIR and 
optical fluxes are corrected for dust extinction using the inter- 
stellar reddening curves of Schild (1977), Cardelli, Clayton, & 
Mathis (1989), and Mathis (1990). A 10% uncertainty is as- 
sumed for the Galactic extinction. 

4.1. GRB 980703 

The emission of the afterglow of GRB 980703 is dominated 
by the host galaxy at only few days after the main event. No 
break has been detected in the optical within this time inter- 
val, which means that we can set only a lower limit on the jet 
opening angle, by requiring tj to be sufficiently large. 

The decay of the optical emission is characterized by a 
power-law index a Q = 1.17 ± 0.25 (Bloom et al. 1998), or 
a a = 1.39 ± 0.30 (Castro-Tirado et al. 1999a), or a a = 
1.63 ± 0.12 (Vreeswijk et al. 1999), thus the average index 
is TxZ — 1.53 ± 0.10. The slope of the optical spectrum is 
Pa = 2.71 ±0.12 (Vreeswijk etal. 1999) at t = 1.2 day, and an 
X-ray decay index a x — 1.33 ± 0.25 was found by Galama et 
al. (1998). Thus the observations give a x -oT = -0.20±0.27. 
In view of the analytical considerations given in §3. 1 , the above 
q x — marginally rules out only a model with s = and 
v c < v x . 

For s = and v a < v c , the optical decay index is a a ~ 
(3p — 3)/4, thus observations require that p = 3.15 ± 0.16, 
which implies (3 a = (p - l)/2 = 1.07 ± 0.08. Such a spec- 
trum is much harder than observed, therefore consistency be- 
tween observations and the fireball model can be achieved only 
if there is a substantial extinction. Given that the Galactic ex- 
tinction toward this afterglow is E(B — V) — 0.061 (Bloom et 
al. 1998), it follows that most this extinction is due to the host 
galaxy. Vreeswijk et al. (1999) have shown that the synchrotron 
power-law spectrum becomes consistent with the observed one 
for an intrinsic extinction of Ay = 1.45 ± 0.13. At the same 



time the dereddened optical fluxes lead to (3 ox = 1.06 ± 0.04, 
which is consistent with the dereddened j3 Q = (p — l)/2, im- 
plying that the cooling break is above X-rays. 

For s = 2 and assuming v Q < v c < v x , leads to a Q = 
(3p — l)/4, therefore the observed a requires p — 2.48 ± 0. 16 
and/3 D = (p— l)/2 = 0. 74 ±0.08, which is again much harder 
than observed. Using the approximations found by Cardelli, 
Clayton, & Mathis (1989) for the UV interstellar reddening, 
one can show that dust extinction in the host galaxy steep- 
ens the synchrotron spectrum by Af3 ~ 1.1 3 Ay at the fre- 
quency which is red-shifted in the observer l/-band. There- 
fore the synchrotron and observed spectra are consistent if 
Ay = 1.74 ± 0.13, which is larger than for the s = model 
because a wind-like medium yields softer spectra for the same 
decay index a. The corresponding dereddened optical emission 
has (3 ox = 1.11 ± 0.04, which falls between the limiting values, 
(p — 1 ) /2 and p/2, allowed by the synchrotron model. 

Thus the afterglow of GRB 980703 may be explained by 
models with homogeneous external media, v a < v x < v c , p ~ 
3.1, and an intrinsic extinction Ay = 1.45 ± 0.13, or by wind 
models with v Q < v c < v x , p ~ 2.5, and Ay — 1.74 ± 0.13. 
Figure 1 shows an s = model yielding an acceptable fit to the 
data. 

4.2. GRB 990123 

After the subtraction of the host galaxy, the r-band light- 
curve of the afterglow of GRB 990123 had a break around 
tb = 2 day (Kulkarni et al. 1999a), with the asymptotic log- 
arithmic slope changing by Aa ~ 0.55 ± 0.07 from that at 
early times, a 0j i = 1.10 ± 0.03, to a ,2 = 1-65 ± 0.06 after tb- 
The same break magnitude is implied by the power-law indices 
found by Castro-Tirado et al. (1999b). The reported slopes of 
the optical spectrum are (3 Q — 0.75 ± 0.23 at t = 1.2 day 
(Galama et al. 1999), and (3 = 0.8±0.1 at t = 1 day (Kulkarni 
et al. 1999a). The optical-to-X-ray spectral slope (3 ox reported 
by Galama et al. (1999) is (3 ox = 0.67 ± 0.02 at t = 1.2 day. 
Therefore (3 ox — f3 a = —0.1 ± 0.1, indicating that the cooling 
break is not between optical and X-ray frequencies at t < 
The first BeppoSAX measurement (Heise et al. 1999) and the 
ASCA data (Murakami et al. 1999) give an X-ray decay index 
a x = 1.17 ± 0.10 for 0.2 d < t < 2d, therefore a x ~ a a ,i, 
consistent with both v c < v Q and v x < v c . 

If the break in the optical emission were due to the passage 
of the v c , then Aa > 1/4 requires s = 2 and an increas- 
ing v c (see §3.1). Equation ( p6| ) and the observed Aa give 
p = 2.46 ± 0.08. Then at t < tt, when v c < v a , the optical- 
to-X-ray slope should be f3 ox = p/2 — 1.23 ± 0.04, clearly 
inconsistent with the observations. This shows that the optical 
break seen in this afterglow is not caused by the passage of v c . 

A cooling frequency below optical cannot be accommo- 
dated by a jet model either, as the observed (3 ox would im- 
ply p = 2f3 ox = 1.34 ± 0.04, leading to a decay index 
c*o,i = (3p — 2)/4 = 0.51 ± 0.03 (irrespective of the type of 
external medium) inconsistent with the observed value. There- 
fore the cooling frequency must be above X-ray, implying 
p = 2j3 ox + 1 = 2.34 ± 0.04. Then s = leads to 
a o l = (3p - 3)/4 = 1.00 ± 0.03, while s = 2 yields 
Oio.i — (3p — l)/4 = 1.50 ± 0.03. The latter case is clearly 
inconsistent with the observations. 

From this analysis one can conclude that a successful model 
for the afterglow of GRB 990123 is that of a jet interacting 
with a homogeneous external medium, with parameters such 
that v D < v x < v c at t < t{,, and p ~ 2.3 . This value of p im- 
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plies an optical decay index a .2 ~ p substantially larger than 
found by Kulkarni et al. (1999a). However, the observations 
made after tj ~ 2 day do not span a sufficiently long time and 
may underestimate the asymptotic a D .2- 

The radio and optical data at t = 1.2 day give a radio-to- 
optical slope f3 ro = 0.27 ± 0.04. If the injection break were 
below radio frequencies (v r ) at this time, then p = 2f3 ro + 1 = 
1.54 ± 0.08, inconsistent with the value derived above from the 
optical and X-ray data. Therefore v r < v. L at t < tj . But in 
this case, as pointed out by Kulkarni et al. (1999a), the radio 
emission should rise until Vi falls below v r . For a jet, this rise 
stops around tj and is followed by a constant emission until 
Vi = v r , after which the radio flux should decrease. How- 
ever, the radio emission of 990123 exhibits a strong dimming 
around 2 days (see Figure 2) which, as suggested by the argu- 
ment above, cannot be explained by the forward shock emis- 
sion in a jet model (see the radio light-curve shown in Figure 
2 with dashed line). We shall assume that the two earliest ra- 
dio fluxes are produced by another radiation mechanism, for 
instance emission by cooled electrons that were accelerated by 
the reverse shock (Sari & Piran 1999) or emission from less rel- 
ativistic ejecta surrounding the jet, and use these fluxes only as 
upper limits for the forward shock emission. 

Figure 2 shows a jet model for the emission of 990123. The 
X-band fluxes observed after 10 days lie above the model pre- 
diction, being inconsistent with an achromatic break resulting 
from collimation of ejecta. If the optical flash of GRB 990123 
was due to a reverse shock propagating in the ejecta (Sari & 
Piran 1999), and if the peak of this flash, seen at t ~ 50 sec- 
onds, corresponds to the fireball deceleration timescale, then 
the isotropic equivalent energy and external density of the 
model in Figure 2 imply that the fireball initial Lorentz fac- 
tor was T = 1400 ± 700. Alternatively, the optical flash 
may have been caused by internal shocks in an unstable wind 
(Meszaros & Rees 1999, Kumar & Piran 2000), in which case 
Tq is more uncertain. 

4.3. GRB 990510 

The optical afterglow of GRB 990510 exhibited a break 
around tb = 1.5 day, whose reported magnitude is Aa = 
1.80 ± 0.20 (Israel et al. 1999) in the ^-band, Aa a = 1.67 ± 
0.02 (Stanek et al. 1999), to Aa a = 1.36 ± 0.05 (Harrison 
et al. 1999). The optical asymptotic decay indices found by 
Harrison et al. (1999) are a ,i = 0.82 ± 0.02 and a a . 2 = 
2.18 ± 0.05. The X-ray decay index was a x = 1.42 ± 0.07 
(Kuulkers et al. 1999) at 0.3d < t < 2d, while the optical 
spectral slope was O = 0.61 ± 0.12 (Stanek et al. 1999) at 
t = 0.9 day. The available data imply an optical-to-X-ray slope 
p ox = 0.90 ± 0.04 at t = 0.72 day. 

Even the smallest reported Aa is above 5 /4, therefore the 
break seen in t he optical emission cannot be due to the passage 
of v c (see §3.3) and must have been caused by jet effects. From 
Pox — Po — 0.29 ± 0.13 one can infer that v c is between optical 
and X-rays. The same conclusion is suggested by a ,i < a x , 
however the X-ray observations were made at times close to 
tj, thus the observed X-ray decay index may overestimate the 
asymptotic a x at earlier times. 

For a homogeneous medium a a .i = {'ip — 3)/4, therefore 
observations require thatp = 2.09 ± 0.03, thus the analytically 
expected values of a .2 ~ P and (3 a = (p— l)/2 = 0.55 ± 0.02 
are consistent with those observed. For a wind-like medium 
&o,i = (3p— l)/4, requiring that p = 1.43 ± 0.03, which leads 
to an optical slope f3 a = (p — l)/2 = 0.22 ± 0.02 inconsis- 



tent even with the softest optical spectrum reported by Stanek 
et al. (1999): j3 = 0.46 ± 0.08. Such a small value of p also 
implies j3 ox < p/2 = 0.72 ± 0.02, again inconsistent with the 
observations. Thus a wind-like medium is ruled out. 

Therefore the afterglow of GRB 990510 can be accommo- 
dated by a model with s = 0, v Q < v c < v x and p ~ 2.1. Note 
that the model shown in Figure 3 fits well the X-ray data, the 
light-curve steepening being very slow. 

4.4. GRB 991216 

The optical decay index of the afterglow of GRB 991216 was 
a 0) i = 1.22 ± 0.04 (Halpern et al. 2000, Sagar et al. 2000) at 
t < t), = 2 day. Halpern et al. have shown that a broken power- 
law optical light-curve with a Q ^ = 1-53 ± 0.05 at t > tb plus 
the contribution R — 24.8 ±0.1 from a galaxy located at ~ 1 
arcsecond from the optical transient explains well the observa- 
tions. Sagar et al. (2000) find that two measurements made 
after tb are dimmer by 2a than the power-law extrapolation 
from earlier fluxes. Therefore there is evidence that the decay 
of the optical emission of the 991216 afterglow steepened by 
Aa = 0.31 ± 0.06 

The X-ray data span 1.3 decades in time before tb and have 
a x = 1.62 ± 0.07 (Frail et al. 2000, Halpern et al. 2000), 
leading to a. x — a .i — 0.40 ± 0.08. The optical spectrum 
is puzzling, exhibiting a turn-over at the J-band at t ~ 1.5 day 
(Frail et al. 2000, Halpern et al. 2000), although the J and K- 
band measurements reported by Garnavich et al. (2000) restore 
a power-law spectrum of slope fi a = 0.58 ± 0.08 at t = 1.7 
day. Garnavich et al. (2000) have pointed out that, if extinction 
is overestimated by a factor 1.3-1.5 close to the Galactic plane 
(Stanek et al. 1999), the dereddened optical spectrum could be 
softer: p o = 0.87±0.08. The optical-to-X-ray spectral slope is 
Pox = 0.80 ± 0.10 (Garnavich et al. 2000, Halpern et al. 2000) 
ati = 1.7 day. 

The faster decay seen in the X-rays than in the optical before 
tb requires v Q < v c < v x and a homogeneous external medium 
(see §[3.3|), in which case it is expected that a x — a a \ < 1/4. 
This difference is 2a below the observed value. For v Q < v c , 
the analytical optical decay index is a 0j i = (3p — 3)/4, thus 
observations imply p = 2.63 ± 0.05. This leads to p o = 
(p — l)/2 = 0.82 ± 0.03, which is consistent with the softer 
dereddened suggested by Garnavich et al. (2000). However, for 
such a spectrum, p ox — Po = —0.07 ± 0.13 does not support a 
cooling frequency between optical and X-rays. 

If the optical light-curve break were due to the passage of 
v c then Aa = 1/4, consistent with the observations. Then at 
t ~ tb, when v c ~ v Q , the optical-to-X-ray slope should be 
Pox — p/2 = 1.32 ± 0.03, clearly inconsistent with the obser- 
vations. Therefore the break in the optical emission of 991216 
was not caused by the passage of v c . Instead, the small magni- 
tude optical break could be explained by jet effects if the short 
time baseline of the observations after tj capture only a part of 
the full steepening. 

However there are some major difficulties that a jet model 
with v Q < v c < v x encounters. The decay index of the radio 
emission a r = 0.82 ± 0.02 (Frail et al. 2000) at t > 1 day re- 
quires the injection break Vi to be below 10 GHz, which leads 
to a radio-to-optical spectral slope p ro = Po — 0.6 -j- 0.9, while 
observations give p ro = 0.20 ± 0.05 at 2 days. Furthermore, 
the decay indices at radio and optical frequencies should be the 
same, yet observations show that a 0l \ — a r = 0.40 ± 0.04. 
The discrepancy is large enough to suggest that a model with 
the above features cannot accommodate the optical and all the 
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radio data. A possible solution is to "decouple" a r and a .i 
by assuming that the radio emission until several days has a 
different origin, as we assumed for the early radio emission of 
990123. Then the quasi-flat behavior of the radio data at 8-50 
days (Figure 4) is suggestive of a jet undergoing lateral spread- 
ing and with Vi above radio until 50 days, as this is the only 
way of obtaining a constant flux if the external medium is ho- 
mogeneous (the analytical prediction for v < vt and after the 
jet edge is seen, F v oc t -1 / 3 , is based on some approxima- 
tions that are not sufficiently accurate; numerically we obtain a 
slightly different behavior F v ~ const). 

There is, however, a model which does not require an- 
other emission mechanism besides the forward shock or a two- 
component structure of the jet (Frail et al. 2000), and which 
can accommodate the entire radio data. The almost constant 
flux exhibited by the radio emission at 1-3 days (Figure 4) can 
be explained by jet expanding laterally before 1 day, with the 
steepening of the radio emission at several days being due to 
the i/;-passage. Then the steepening of the optical decay at 
few days is not due to jet effects, as suggested by Halpern et 
al. (2000), but to the passage of a spectral break. The cool- 
ing frequency v c does not offer a self-consistent picture: it 
should be above optical shortly before 2 days, as required by 
a x > a .i and, according to equation (28), it should increase 
in time because during the jet sideways expansion phase, thus 
it cannot cross the optical domain. 

Within our afterglow modeling there is only one remaining 
possibility: the steepening of the optical emission is due to 
the passage of the I'm frequency associated with the high en- 
ergy break of the electron distribution (eq. JlO|]). As the jet 
model yields F v oc t~ p at all frequencies above z/j (eq. p^]), 
the observed a 0) i implies a very hard electron distribution with 
p ~ 1.2. Then the observed (3 — 0.58 ± 0.08 requires v c < v a 
at 2 days, to obtain consistency with the analytical expectation 
/3 o = P /2~0.6. 

Such a jet model with a hard electron distribution and cool- 
ing frequency below optical is shown in Figure 4. The observed 
optical emission steepens only mildly after 2 days, therefore 
the high energy break at jm is not too strong. For simplicity 
we have approximated this break as a softening of the electron 
distribution from 7~ p to ^~(p+ s p) at 7m, with Sp a free param- 
eter. The location of v m is set by the fractional electron energy 
e = 0.5 up to 7m (eq- JToj]). Smaller values of e, but larger 
than 0.1, also provide acceptable fits. The afterglow shown in 
Figure 1 is mildly relativistic after 10 days, so that departures 
from the analytically expected light-curve F v>Vi oc t~ p allow 
the model to accommodate both the t~ 12 optical decay before 
2 days and the t~ 8 fall-off of the radio emission after 10 days. 
As illustrated in Figure 4 the effect of interstellar scintillation 
(Goodman 1997) is essential in explaining the departures be- 
tween the observations and the model radio fluxes. 

We note that a jet interacting with a wind-like external 
medium with A* ~ 1 yields fits with acceptable % 2 , but pro- 
duces millimeter fluxes that exceed the 2a upper limits shown 
in Figure 4 (right panel). 

4.5. Parameter Ranges and Afterglow Energetics 

The electron index p and the initial jet opening angle 8q are 
determined by the index of the power-law emission decay and 
by the time when jet effects set in, respectively (see eqs. [|27|l, 
[[2l|], and [^2||). The remaining four model parameters, £q, n, 
e e , and sb, can be determined from the three break frequen- 
cies v a , vi, v c (see eqs. [§], Jl3|], JT^], and [|l7|]), and the syn- 



chrotron flux at the peak of the spectrum. Finding the location 
of the spectral breaks requires observations spanning a wide 
frequency range, from below the lowest break (self-absorption) 
to above the highest break (cooling, more likely). Even if this 
requirement is satisfied, observations in only three domains (ra- 
dio, near infrared-optical, and X-ray) do not determine all the 
spectral breaks unless the unconstrained break(s) cross(es) an 
observing frequency. This does not seem to be the case for the 
afterglows analyzed here, observations providing at most three 
"strong" constraints for four free model parameters. The num- 
ber of constraints is even smaller if the interval between two 
adjacent observational frequencies does not contain a spectral 
break. For instance, in the case where v x < v c , the AT -ray 
fluxes can be predicted from the optical emission and the spec- 
tral slope 0{p). If consistency is found between the model and 
the observed A-ray fluxes, i.e. if the condition v x < v c is in- 
deed satisfied, then the A-ray data will provide only a "weak" 
constraint, as they set only a lower bound on the cooling break 
frequency. 

Therefore the uncertainty in the parameters of the afterglows 
analyzed in this work arises from that the number of effective 
observational constraints is smaller than the number of model 
parameters. To assess these uncertainties, we find sets of pa- 
rameters which yield acceptable fits to the data, for a given 
isotropic equivalent energy £q. The distributions of the parame- 
ters for which the probability of exceeding the % 2 of the data is 
at least 20% are shown in Figures 5 and 6. Note that 8q ranges 
from about 1° to 4° and that the external medium density spans 
four decades. Figure 6 shows the lack of a pattern in the elec- 
tron and magnetic field parameters among the four afterglows. 
As expected, the index p is well constrained for each afterglow, 
the light-curve decay being very sensitive to it. Note that, for 
the afterglows modeled here, p does not have a universal value. 

Figure 7 shows the distributions of jet energy Eq = 
(£q/2)(1 — cos^o) and of the burst 7-ray efficiency, defined 
as the ratio of the isotropic equivalent £~ t of the energy released 
by the GRB in the 20 keV - 1 Me V range to the total fireball en- 
ergy £ 1 + £0. The £ 1 was calculated from the reported fluences 
<£> 7 above 20 keV (Kippen et al. 1998-9) and the measured red- 
shifts: £ 1 ~ 10 53 erg for GRB 980703, £ 7 - 3 x 10 54 erg 
for GRB 990123, £ 7 ~ 2 X 10 53 erg for GRB 990510, and 
£ 7 - 7 x 10 53 erg for GRB 991216. Note that, when O can be 
determined from observations, the jet energies are clustered in 
the 10 50 to 10 51 erg range, and that the minimum efficiency for 
GRBs 990123 and 991216 exceeds 20%. 



5. CONCLUSIONS 

We model the emission of GRB afterglows in the framework 
of collimated, uniform, lateral spreading jets interacting with 
an external medium. The model was used to determine the ini- 
tial jet energy Eq, opening angle Qq, external medium density 
n, and parameters e e and e b quantizing the minimum random 
Lorentz factor of shock-accelerated electrons and the strength 
of the magnetic field, respectively, for the afterglows of GRB 
980703, 990123, 990510, and 991216. 

As illustrated in Figures 5 and 6, the jet aperture 8q and the 
index p of the power-law distribution of electrons injected in the 
down-stream region are well constrained by observations, as the 
effects of collimation are seen at a time which depends strongly 
P), while p determines the afterglow decay rate 
J). Other model parameters - Eq, n, and sb 
- are less well determined as observations in three frequency 
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domains (radio, optical, and X-ray) provide at most three con- 
straints. Observations in a fourth domain, sub-millimeter, mil- 
limeter, or far-infrared, could help determine all the afterglow 
parameters, provided that the spectral breaks are located be- 
tween adjacent observational frequencies. 

Analysis of the available data for the 990123 and 990510 af- 
terglows and the jet interpretation of the break exhibited by their 
optical emissions rule out a wind-like profile for the medium 
which decelerates the jet. However, the emission of the 980703 
and 991216 afterglows can be accommodated by both types 
of external media (homogeneous or a pre-ejected wind). For 
980703, 990123, and 990510, the allowed range of external 
densities is below 1 cm~ 3 , indicating that these bursts did not 
occur in hydrogen clouds. 

For those afterglows with optical light-curve breaks we find 
jet energies lying mostly in the 10 50 -10 51 erg range (Figure 7), 
and jet initial half-angles below 4°. For the afterglow of GRB 
970508, the only one observed in radio, optical and X-ray that 
is not included in this work due to its non-standard optical and 
X-ray brightening, Frail, Waxman & Kulkarni (2000) found a 



similar jet energy, Eq = 5 x 10 erg, a much larger initial 
half-angle, Oq ~ 30°, and an external density n ~ 1 cm -3 . 

The minimum BATSE range (20 keV - 1 MeV) efficiency 
of the afterglows whose optical behavior indicates that the 
ejecta was well collimated ranges from 3% to 30% (Figure 
7). The former limit is within the reach of current calculations 
of internal shock efficiency (Kumar 1999, Lazzati et al. 1999, 
Panaitescu, Spada & Meszaros 1999), but the latter exceeds it. 
However, if the energy distribution within the jet aperture is far 
from isotropy, such that the GRBs we see have an energy-per- 
solid angle peaking in the direction toward the observer, then 
the minimum required efficiency can be significantly smaller. 
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FIG. 1 . — Light-curves (left panel) and t = 1.20 day spectrum (right panel) for the afterglow of GRB 980703 (z = 0.966). Radio data are from Frail et al. (1999). 
We assume 20-40% uncertainty when no errors were reported, the larger errors being for earlier or lower frequency observations, the smaller for later or higher 
frequencies. Optical fluxes are calculated from the magnitudes published by Bloom et al. (1998), Castro-Tirado et al. (1999a), and Vreeswijk et al. (1999). The host 
galaxy contribution (Vreeswijk et al. 1999) is subtracted. Optical fluxes are corrected for a Galactic extinction of E(B — V) = 0.061 (Bloom et al. 1998) and for 
the host galaxy extinction inferred by Vreeswijk et al. (1999): Ay = 1.45 ± 0.13. X-ray fluxes are calculated from the 2-10 keV band fluxes reported by Vreeswijk 
et al. (1999). Triangles indicate 2<r upper limits. Right panel: optical measurements closest to t = 1.20 day are extrapolated to this time using F v oc t~ a , with the 
indices a reported by Vreeswijk et al. (1999). Open symbols for radio fluxes are extrapolations from the 4.1 and 5.3 days measurements of Frail et al. (1999) with 
the analytically expected behavior F v oc Z 1 / 2 . The model shown has \ 2 = 20 for 21 df and the following parameters: £n = 2.9 X 10 54 erg (isotropic equivalent 
energy), n = 7.8 X 10~ 4 cm -3 , e e = 0.075, sb = 4.6 X 10 — 4 , p = 3.08 . The initial jet angle is lower bounded by the condition that no effects of collimation 
are seen until the later available data, which leads to 9q 2.7°. Radiative losses amount to 1.6% at 10 days. Synchrotron emission is the main electron radiative 
cooling mechanism at all times shown here. 
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FIG. 2. — Light-curves and t = 1.24 day spectrum for the afterglow of GRB 990123 (z = 1.61). Radio data are from Kulkarni et al. (1999b). Optical fluxes 
are calculated from the magnitudes reported by Kulkarni et al. (1999a) by subtracting the host galaxy contribution F r = 0.5 ± 0.1 fijy and K = 22.1 ± 0.3, and 
from the HST measurement reported by Odewahn et al. (1999). Fluxes are corrected for a Galactic extinction of A r = 0.04 (Kulkarni et al. 1999a). X-ray fluxes are 
calculated from the 2-10 keV fluxes reported by Heise et al. (1999) and Murakami et al. (1999), using the quoted spectral indices and assuming a 10% error. Other 
data, used for the spectrum shown in the right panel, are from Galama et al. (1999) and Castro-Tirado et al. (1999b), and are extrapolated to t = 1.24 day using the 
optical decay indices given in Castro-Tirado et al. (1999b). Triangles indicate 2<r upper limits. The model shown has x 2 = 16 for 16 df, excluding the first two radio 
measurements and the group of late A'-band observations (which would increase \ 2 by A\ 2 = 36 for), and parameters: Eq = 4.6 X 10 50 erg (initial jet energy), 
do = 2.0°, n = 4.7 X 10~ 4 cm -3 , e e = 0.13, eb = 1-3 X 10~ 4 , p = 2.32 . The steepening of the optical decay is due to jet effects. About 9% of the afterglow 
energy is radiated until 30 days. Inverse Compton scatterings dominate the electron radiative cooling until several days. The model shown with dashed line yields a 
good fit to the optical and X-ray data, but matches only the first two radio measurements, over-predicting the late time radio emission. 
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FIG. 3. — Light-curves and t = 0.72 day spectrum for the afterglow of GRB 990510 (z = 1.62). Radio data are from Harrison et al. (1999). The triangle 
indicates a 2a upper limit. Optical fluxes are obtained from the magnitudes reported by Beuermann et al. (1999), Fruchter et al. (1999), Harrison et al. (1999), 
Israel et al. (1999), Pietrzynski & Udalski (1999), and Stanek et al. (1999). These fluxes are corrected for dust extinction with E(B — V) = 0.20 (Harrison et 
al. 1999, Stanek et al. 1999). The uncertain, small contribution of the host galaxy, V g <; 28.0 (Bloom et al. 2000, Fruchter et al. 2000), is ignored. X-ray fluxes 
are calculated from the 2-10 keV band fluxes and spectral slope published by Kuulkers et al. (2000). The model shown has x 2 = 42 for 69 df, and parameters: 
Eo = 3.0 X 10 50 erg, 6»o = 2.7°, n = 0.14 cm -3 , e e = 0.046, e B = 8.6 X 10~ 4 , p = 2.01. The steepening of the optical decay is due to the collimation of 
ejecta. The radiative losses amount to 53% at 40 days. At all times shown here, inverse Compton scatterings is the dominant electron radiative cooling mechanism. 
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FIG. 4.— Light-curves and t = 1.52 day spectrum for the afterglow of GRB 991216 (z = 1.02). Radio data are from Rol et al. (1999) and Frail et al. (2000). 
Optical fluxes are calculated from the magnitudes published by Garnavich et al. (2000), Halpern et al. (2000), Sagar et al. (2000), and Schaefer (2000). We subtracted 
the contamination from a galaxy located near the OT, estimated by Halpern et al. (2000) at R = 24.8 ±0.1 (the host galaxy, identified by Vreeswijk et al. 2000, has 
R = 26. 9 ±0.4). Optical fluxes are corrected for a Galactic extinction of E(B — V) = 0.63 (Garnavich et al. 2000, Halpern et al. 2000). X-ray fluxes are calculated 
from the 2—10 keV fluxes reported by Corbet et al. (1999), Piro et al. (1999), and Takeshima et al. (1999). Right panel: the measurements closest to t = 1.52 day 
are extrapolated to this time using the power-law scaling found by Garnavich et al. (2000) in the near-infrared and by Halpern et al. (2000) in the optical. Inset: J 
and K measurements of Halpern et al. (2000) are shown with open symbols, while those reported by Garnavich et al. (2000), consistent with a power-law spectrum, 
are indicated with filled symbols. Triangles indicate 2<r upper limits. The model shown has \ 2 = 37 for 43 df and parameters Eo = 0.89 X 10 50 erg, 0o = 3.4°, 
n = 1.6 cm -3 , e e = 0.024, eb = 0.073, p = 1.36, Sp = 0.7. The vm frequency is set by the fractional electron energy, assumed e = 0.5. The steepening 
of the X-ray and optical fall-off is due to the passage of the vm break. Vertical dotted lines indicate the amplitude of the interstellar scintillation, which we model 
following the treatment given by Walker (1998). The afterglow radiates 82% of its energy until 80 days. The electrons cool radiatively mainly through synchrotron 
emission after few days. 



11 




12 3 4 

6„ (deg) 



980703 


I I 


1 1 


990123 






990510 






" 991216 




1 | 


- 


1 


1 

1 




j 

1 


1 




1 


' !' ^ : 


r'r', i i i i > 


TT^ , l 


i 1 [ i i 1 i r 


-4 


-2 


2 



log n (cm 3 ) 



FIG. 5. — Distributions of the parameters 8q (jet initial opening) and n (external medium density) obtained from fits to the afterglows of GRB 980703, 990123, 
990510, and 991216. For the models with these parameters the probability of exceeding the \ 2 of the data is at least 20%. The decay of the afterglow of GRB 980203 
does not show evidence for collimation of ejecta until several days. The initial jet apertures shown in this case are obtained by requiring that tj = 5 day (eq. ||27|]), 
therefore they are lower limits of the true 6o . 




FIG. 6. — Distributions of the parameters e e (for the minimum energy of the injected electrons), eb (fractional energy in magnetic field), and p (exponent of the 
injected power-law distribution) for the same models as in Figure 5. Note that these four afterglows have distinct allowed ranges for p. 




FIG. 7. — Distribution of the jet energy Eq and GRB radiative efficiency for the models whose parameters are shown in Figures 5 and 6. The radiative efficiency 
is the ratio of the energy £ 7 released by the GRB in the BATSE range (approximately 20 keV - 1 MeV), assuming spherical symmetry, and the initial fireball energy 
£fl, = S 1 + So, where So is the isotropic equivalent of the ejecta energy which yields acceptable fits to the afterglow emission. For GRB 980703 we used the lower 
limits on 8q shown in Figure 5 to calculate the jet energy. 



